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(57) ABSTRACT 


Provided are methods of using electromagnetic waves for 
detecting metal and/or dielectric objects. Methods include 
directing microwave and/or mm wave radiation in a prede- 
termined direction using a transmission apparatus, including 
a transmission element; receiving radiation from an entity 
resulting from the transmitted radiation using a detection 
apparatus; and generating one or more detection signals in 
the frequency domain using the detection apparatus. Meth- 
ods may include operating a controller, wherein operating 
the controller includes causing the transmitted radiation to 
be swept over a predetermined range of frequencies, per- 
forming a transform operation on the detection signal(s) to 
generate one or more transformed signals in the time 
domain, and determining, from one or more features of the 
transformed signal, one or more dimensions of a metallic or 
dielectric object upon which the transmitted radiation is 
incident. A system and method for remote detection and/or 
identification of a metallic threat object using late time 
response (LTR) signals is also disclosed. 
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REMOTE DETECTION AND 
MEASUREMENT OF OBJECTS 


TECHNICAL FIELD 


[0001] The present invention relates to the detection of 
objects, and more particularly, to techniques for remote 
detection and measurement of objects. 


BACKGROUND ART 


[0002] Itis well known to use electromagnetic radiation to 
detect the presence of objects (e.g. handheld detectors used 
for detecting objects on or under the ground, and walk- 
through arches at airports). 

[0003] However, the conventional detectors used at air- 
ports maybe unable to determine the dimensions of objects 
to any significant degree, and thus may be unable to distin- 
guish between objects of different types, i.e. harmless (belt 
buckles, cameras), and potentially dangerous (guns, knives). 
[0004] The detection of concealed weapons, especially 
handguns, may be a very great problem for security appli- 
cations that currently cannot be policed without a non- 
portable system, for example random checks in an urban 
environment. The use of microwaves (electromagnetic 
waves with wavelengths in the centimeter to millimeter 
range) may provide a means for the standoff detection and 
identification of concealed conducting items such as hand- 
guns and knives, Large metal objects, such as handguns, 
may give a significantly different and generally larger 
response when irradiated by low power microwaves than 
that from the human body, clothing and/or benign normally- 
carried objects. The larger response may be detected using 
a combination of antenna and sensitive receiver. 

[0005] By actively illuminating, an object with wide-range 
swept and/or stepped frequency microwave and/or millime- 
ter wave radiation, the frequency response of the return 
signal may give the range and/or information regarding 
dimensions of the object. This method may be substantially 
equivalent to using a fast microwave pulse and measuring 
the response as function of time, as used in conventional 
RADAR. Selecting a part of the return signal within a 
particular range may aid the positive identification of the 
suspect object and may also help to reject background 
signals. The analysis of the time response may give further 
information as to the dimensions of the target. This tech- 
nique may also be applied to the detection of dielectric 
layers, such as, for example, an explosive vest strapped to a 
suicide bomber (see Active millimeter wave detection of 
concealed layers of dielectric material, Bowring N. J., 
Baker J, G., Rezgui M., Southgate M., Proceedings of the 
SPIE 6540-52 2007; and A sensor for the detection and 
measurement of thin dielectric layers using reflection of 
frequency scanned millimetric waves, Bowring N. J., Baker 
J. G., Rezgui R, Alder J. F. Meas. Sci. Technol. 19 024004 
(7pp) 2008). However, such techniques have not been here- 
tofore used for detecting and measuring metal objects. 
[0006] A system based on swept frequency RADAR has 
been proposed (U.S. Pat. Nos. 6,359,582 and 6,856,271 and 
published application US2007/0052576), In the disclosed 
systems, the frequency may be swept by typically by 1 GHz 
around about 6 GHz. The depth resolution that is achievable 
is therefore only 15 cm, thus the system may not give details 
of the objects, The detection relies on comparing gross 
features of the signal as a whole with similar suspicious and 


Dec. 26, 2019 


benign signals to which the system had been previously 
exposed, Also the measurement of polarization properties of 
the scattered signal may be used. 

[0007] In the aforementioned patent documents, the low 
frequency of operation makes the angular resolution of the 
antennae poor and the wide field of view makes it difficult 
to single out particular targets and/or to determine on which 
part of the target the threat is situated, This may be improved 
by changing to higher frequencies where microwave optics 
becomes effective. This may be particularly important for 
explosives detection where the contrast from the body signal 
is low, Systems working at higher frequencies but still with 
a limited bandwidth have been proposed by Gorman et al 
(U.S. Pat. No. 6,967,612) and by Millitech (U.S. Pat. No. 
5,227,800). Many systems have been produced to enable 
images of the target to be obtained using either active 
microwave illumination or the passive thermal emission of 
the target (SPIE 2007). These systems use multi-detector 
arrays and some form of mechanical scanning. Passive 
systems, though giving more realistic images, tend to be 
slow and show poor contrast for dielectric targets. Active 
illumination systems can be acquired faster, but may suffer 
from strong reflections from benign objects such as the 
human body, which make it difficult to distinguish from 
metal threat objects, All scanning systems may require 
complex human or Artificial Intelligence interaction to inter- 
pret the image and/or to pick out the suspect features. This 
makes their deployment in many applications difficult. 
[0008] Itis apparent that systems which can identify threat 
objects at standoff distances may have many applications, 
where conventional metal detector booths are inappropriate. 
These may include covert surveillance and mobile operation 
in streets and buildings. Portable, compact and cost-effective 
systems are not presently available and this invention seeks 
to address this need. 


SUMMARY OF THE INVENTION 


[0009] According to some embodiments of the present 
invention there are provided systems for remote detection of 
one or more dimensions of a metallic or dielectric object. 
Some embodiments of such systems may include a trans- 
mission apparatus, including a transmission element, for 
directing microwave and/or mm wave radiation in a prede- 
termined direction, a detection apparatus, for receiving 
radiation from an entity resulting from the transmitted 
radiation and generating one or more detection signals in the 
frequency domain, and a controller. In some embodiments, 
the controller may be operable to (1) cause the transmitted 
radiation to be swept over a predetermined range of fre- 
quencies, (11) perform a transform operation on the detection 
signal(s) to generate one or more transformed signals in the 
time domain or optical depth domain, and (111) determine, 
from one or more features of the transformed signal, one or 
more dimensions of a metallic or dielectric object upon 
which the transmitted radiation is incident. 

[0010] In some embodiments, the transmission element is 
a directional element that may be pointed by a user. 
[0011] In some embodiments, the controller is operable to 
initiate step (1) upon receiving an activation signal, the 
activation signal corresponding to a user input and/or to 
detection of the presence of the entity. 

[0012] In some embodiments, step (1) comprises stepwise 
sweeping by predetermined steps in frequency and step (11) 
comprises (iia) performing a transform operation-after each 
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sweep to produce a time domain or optical depth domain 
trace, and (iib) storing each time domain trace in a respective 
sweep channel, said time domain or optical depth domain 
traces thereby comprising said transformed signals. 


[0013] In some embodiments, step (111) comprises normal- 
ising said according to range one or more transformed 
signals. 


[0014] In some embodiments, step (iii) further comprises 
using a Complex Fourier Transform and/or Direct Fourier 
Transform to convert transformed signals to the x-dimen- 
sion, and determining the x-positions of peaks on the trans- 
formed signals. 


[0015] In some embodiments, step (iii) further comprises, 
from the x-positions, using 


c 


= TAF 


where L=optical depth 
[0016] c=the speed of light 
[0017] Af=the periodicity in the frequency domain, 


to determine corrected x-axis positions, and thereby optical 
depth. 


[0018] In some embodiments, step (ii) includes the pro- 
cedure of Appendix B.1, thereby producing first and second 
outputs (Outputl, Output2) dependent upon the detection 
signals, wherein Output 1 is the sum of all correlations 
between vectors in a first array, the vectors in the first array 
comprising, for each sweep channel, a stored signal above a 
threshold that are derived by Direct Fourier Transform from 
the detection signals, and Output2 is the sum of integrated 
signals above the threshold for each sweep channel. 


[0019] In some embodiments, step (ii) includes the pro- 
cedure of Appendix B.2, thereby producing a third output 
(Output3) dependent upon the detection signals, wherein 
Output3 1s, for each sweep channel, the best (lowest) cor- 
relation value between the ideal response for a number of 
barrel lengths and for a number of weapon calibers stored in 
memory and the direct untransformed detection signals 
IER. 


[0020] In some embodiments, the detection apparatus 
includes a first detection element directed in a first direction, 
towards the entity, for generating non-polarized detection 
signals, and a second detection element, directed at 90 
degrees to the first element, for generating cross-polarized 
detection signals, and wherein step (11) includes the proce- 
dure of Appendix B.3, thereby producing fourth and fifth 
outputs (Output4, Output5) dependent upon the detection 
signals, wherein Output4 is, for each sweep channel, the sum 
of correlations between the transformed signals, the trans- 
formed signals comprising a Complex Fourier Transform of 
the non-polarized detection signals and of the cross-polar- 
ized detection signals, and Output5, for each sweep channel, 
the sum of the integrated non-polarized and cross-polarized 
signals after Complex Fourier Transform. In some embodi- 
ments, the integrated non-polarized and cross-polarized sig- 
nals comprise integrations of transformed signals within one 
or more distance windows, the contents of each distance 
window being stored in a database in association with a 
respective response for a particular weapon. 
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[0021] In some embodiments, step (ii) includes the pro- 
cedure of Appendix B.4, thereby producing sixth and sev- 
enth outputs (Output6, Output7) dependent upon the detec- 
tion signals. 

[0022] In some embodiments, the system further includes 
a neural network, the neural network having as inputs 
thereto any combination of (a) the first and second outputs 
(Output 1, Output2), (b) the third output (Output3), (c) the 
fourth and fifth outputs (Output4, Outputs), and (d) the sixth 
and seventh outputs (Output6, Output7), wherein the output 
of the neural network is an indication of a confidence level 
of a metallic or dielectric object of a predetermined type 
being detected, for example, 1=gun detected, 0—no metallic 
or dielectric object detected. 

[0023] According to some embodiments of the present 
invention there are provided methods for remote detection of 
one or more dimensions of a metallic or dielectric object. 
Such methods may include a transmission apparatus, includ- 
ing a transmission element, for directing microwave and/or 
mm wave radiation in a predetermined direction, a detection 
apparatus, for receiving radiation from an entity resulting 
from the transmitted radiation and generating one or more 
detection signals in the frequency domain, and a controller, 
Methods may include operating the controller to (1) cause 
the transmitted radiation to be swept over a predetermined 
range of frequencies, (ii) perform a transform operation on 
the detection signal(s) to generate one or more transformed 
signals in the time domain, and (iii) determine, from one or 
more features of the transformed signal, one or more dimen- 
sions of a metallic or dielectric object upon which the 
transmitted radiation is incident. In some embodiments, the 
transmission element is a directional element that may be 
pointed by a user. 

[0024] 11) some embodiments, the controller is operable to 
initiate step (i) upon receiving an activation signal, the 
activation signal corresponding to a user input or to detec- 
tion of the presence of the entity. 

[0025] In some embodiments, step (1) comprises stepwise 
sweeping by predetermined steps in frequency and step (11) 
comprises (iia) perform a transform operation after each 
sweep to produce a time domain or depth domain trace, and 
(iib) storing each time domain or depth domain trace in a 
respective sweep channel, said time domain or depth domain 
traces thereby comprising said transformed signals. 

[0026] In some embodiments, step (111) comprises normal- 
ising according to range one or more transformed signals. 

[0027] In some embodiments, step (111) further comprises 
using a Complex Fourier Transform and/or Direct Fourier 
Transform to convert transformed signals to the x-dimen- 
sion, and determining the x-positions of peaks on the trans- 
formed signals. 

[0028] In some embodiments, step (111) further comprises, 
from the x-positions, using 


where L-distance to entity 
[0029] c=the speed of light 
[0030] Af=the periodicity in the frequency domain, 


to determine corrected x-axis positions, and thereby optical 
depth. 
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[0031] In some embodiments, step (ii) includes the pro- 
cedure of Appendix B.1, thereby producing first and second 
outputs (Outputl, Output2) dependent upon the detection 
signals, wherein Outputl is the sum of all correlations 
between vectors in a first array, the vectors in the first array 
comprising, for each sweep channel, a stored signal above a 
threshold that are derived by Direct Fourier Transform from 
the detection signals, and Output2 is the sum of integrated 
signals above the threshold for each sweep channel. 


[0032] In some embodiments, step (ii) includes the pro- 
cedure of Appendix 13.2, thereby producing a third output 
(Output3) dependent upon the detection signals, wherein 
Outputs is, for each sweep channel, the best (lowest) cor- 
relation value between the ideal response for a number of 
barrel lengths and for a number of weapon calibers stored in 
memory and the direct untransformed detection signals 
IE4I?. 

[0033] In some embodiments, the detection apparatus 
includes a first detection element directed in a first direction, 
towards the entity, for generating non-polarized detection 
signals, and a second detection element, directed at 90 
degrees to the first element, for generating cross-polarized 
detection signals, and wherein step (ii) includes the proce- 
dure of Appendix B.3, thereby producing fourth and fifth 
outputs (Output4, Outputs) dependent upon the detection 
signals, wherein Output4 is, for each sweep channel, the sum 
of correlations between the transformed signals, the trans- 
formed signals comprising a Complex Fourier Transform of 
the non-polarized detection signals and of the cross-polar- 
ized detection signals, and Output5 is, for each sweep 
channel, the sum of the integrated non-polarized and cross- 
polarized signals after Complex Fourier Transform, In some 
embodiments, the integrated non-polarized and cross-polar- 
ized signals comprise integrations of transformed signals 
within one or more distance windows, the contents of each 
distance window being stored in a database in association 
with a respective response for a particular weapon. 


[0034] In some embodiments, step (ii) includes the pro- 
cedure of Appendix B.4, thereby producing sixth and sev- 
enth outputs (Output6, Output7) dependent upon the detec- 
tion signals. 


[0035] In some embodiments, methods may further com- 
prises providing a neural network, the neural network hav- 
ing as inputs thereto any combination of (a) the first and 
second outputs (Outputl, Output2), (b) the third output 
(Output3), (c) the fourth and fifth outputs (Output4, Out- 
puts), and (d) the sixth and seventh outputs (Output6, 
Output7), wherein the output of the neural network is an 
indication of a confidence level of a metallic or dielectric 
object of a predetermined type being detected, for example, 
1=gun detected, 0—no metallic or dielectric object detected. 


[0036] Some embodiments of the present invention may 
include a recorded or recordable medium having recorded or 
stored thereon digital data defining or transformable into 
instructions for execution by processing circuitry, the 
instructions corresponding the systems and methods herein. 


[0037] Some embodiments of the present invention may 
include a server comprising processing circuitry, memory 
and a communications device, the server being programmed 
for communicating on demand or otherwise digital data 
defining or transformable into instructions for execution by 
processing circuitry, the instructions corresponding to the 
systems and methods herein. 
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[0038] Techniques according to some embodiments of the 
invention entail actively illuminating an object with wide 
range stepped microwave or millimeter wave radiation and 
inducing a local electromagnetic field on the surface of the 
object and within the barrel, comprised of a superposition of 
modes. 

[0039] The coupling to these modes from the illuminating 
and scattered fields is, in general, frequency dependent and 
this forms the basis for the detection and identification of 
conducting items. The object may be fully illuminated if a 
full spectrum of modes and therefore a full frequency 
response are to be excited and collected. 

[0040] The scattered EM power may be typically mea- 
sured at “stand off" distance of several meters as the 
illuminating field is frequency swept over as wide a range as 
possible and patterns in frequency response characteristic to 
the target object being sought are looked for. 

[0041] This system may rely on contributions from one or 
more of the following effects. 

(1) Swept reflectrometry, in which the distances between 
corners, edges and cavities on the weapon independently 
from the distances to the source and the detector are ascer- 
tained in real-time. 

(2) Barrel tone detection. The caliber and barrel length ofthe 
weapon can be ascertained when the barrel is oriented 
towards the general direction of a detector and source 
(+45”). This is determined by an aspect invariant method 
outlined in the appropriate section below. 

(3) Cross-polarization effects. The use of two or more 
detectors, the first oriented in the same direction as the 
illuminating radiation and the second at 90? (the cross 
polarized detector) can yield important effects, The time 
dependent responses in the cross polarized detector are 
enhanced when a handgun is present. Furthermore, if the 
frequency is swept over a large range, the reflections from 
various corners of the object can be resolved in the time 
domain, These reflections are different or anti-correlated in 
the cross polarized detector. This differentiates the technique 
from previously known polarization-based detectors out- 
lined above. 

(4) Aspect Independent effects (Late Time Responses). A 
late time response or recurrence from the object is identified, 
to effectively give information related to the size and dimen- 
sions of the object. This enables dimensions to be deter- 
mined in truly aspect independent manner. This phenom- 
enon is the further set of scattered signals obtained after the 
main radiation has returned to the receiver, assuming the 
application of a pulse of radiation, It is caused by locally 
induced currents in the object re-radiating. The timing of the 
response can be directly linked to the size and natural 
resonances of the object, and this is illustrated herein by 
simple objects with one or more dimensions. It can, there- 
fore, prove a useful tool in the identification of weapons, 
optionally when combined with the other techniques iden- 
tified herein, 

[0042] The concept of target identification by utilising 
re-radiation from a metal object illuminated by a radar pulse, 
after that pulse has passed the object, has been used and 
understood for decades, primarily in the field of radar for 
aircraft and missile identification purposes [e.g. Baum C. E 
*On the singularity expansion method for the solution of 
electromagnetic interaction problems", Air Force Weapons 
Lab. Interaction Notes, Note 88, December, 1971] The 
present inventors demonstrate in this disclosure how tech- 
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niques that build upon the basic principles can be used to 
confirm the presence or the absence of weapons, especially 
handguns and knives, concealed on or about the human 
body. 


[0043] When a highly conductive object, such as a metal 
handgun or knife, is illuminated by a sudden pulse of 
microwave frequency electromagnetic radiation there are 
surface currents excited on that object. When the exciting 
pulse has passed the object, the excited surface currents 
oscillate in such a way as to give rise to re-radiation that 
carries an electromagnetic signature which is unique to that 
object. Of critical importance and underpinning the tech- 
niques according to some embodiments is the fact that the re 
radiated signature is independent of such vagaries as pulse 
shape and object (target) orientation. This re-radiated elec- 
tromagnetic signal is known as the Late Time Response 
(LTR) and it is this signal that is excited, captured and 
processed to interrogate the human body for concealed 
weapons, according to some embodiments of the present 
invention, 


[0044] The techniques identified in this application are 
suitable for a deployable gun and concealed weapons detec- 
tion system and does not rely on imaging techniques for 
determining the presence of a gun. One or more of the 
techniques can be combined to reduce the false-positive 
events from the detector. Hereinafter, experimental sets of 
responses from typical metal or partially conducting objects 
such as keys, mobile phones and concealed handguns at a 
range of frequencies are presented, According to another 
aspect of the invention there is provided a system for remote 
detection and/or identification of a metallic threat object, 
comprising: a transmission apparatus, including a transmis- 
sion element, configured to direct microwave and/or mm 
wave radiation in a predetermined direction, a detection 
apparatus configured to receive radiation from an entity 
resulting from the transmitted radiation and to generate one 
or more detection signals in the frequency domain, and 
control circuitry, the control circuitry being operable to (1) 
cause the transmitted radiation to be swept over a predeter- 
mined range of frequencies, (11) perform a transform opera- 
tion on the detection signal(s) to generate one or more 
transformed signals in the time domain, (111) extract, from 
one or more features of the transformed signal, the late time 
response (LTR) signal of the metallic threat object upon 
which the transmitted radiation is incident, the metallic 
threat object being carried by or associated with the entity; 
and (iv) from the LTR signal, determine the presence and/or 
identity of the metallic threat object. According to another 
aspect ofthe invention there is provided a method for remote 
detection and/or identification of a metallic threat object, 
comprising: providing a transmission apparatus, including a 
transmission element, configured to direct microwave and/or 
mm wave radiation in a predetermined direction, providing 
a detection apparatus configured to receive radiation from an 
entity resulting from the transmitted radiation and to gen- 
erate one or more detection signals in the frequency domain, 
and control circuitry, the method comprising (1) causing the 
transmitted radiation to be swept over a predetermined range 
of frequencies, (ii) performing a transform operation on the 
detection signal(s) to generate one or more transformed 
signals in the time domain, (iii) extracting, from one or more 
features of the transformed signal, the late time response 
(LTR) signal of the metallic threat object upon which the 
transmitted radiation is incident, the metallic threat object 
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being carried by or associated with the entity; and (iv) from 
the LTR signal, determine the presence and/or identity of the 
metallic threat object, According to another aspect of the 
invention there is provided a threat object detection system, 
for remote detection and/or identification of a metallic threat 
object being carried by or associated with the entity, com- 
prising: a plurality of detection systems, each detection 
system being in accordance with any of claims 28 to 59 of 
the appended claims, the detection systems being arranged 
spaced apart in relation to a closed path defined by walls and 
to be traversed by the entity, whereby, in use, transmitted 
radiation from at least one of the detection systems is 
incident upon the entity, According to another aspect of the 
invention there is provided a computer program product for 
remote detection of one or more dimensions of a metallic or 
dielectric object, the computer program product comprising 
a computer usable storage medium having computer read- 
able program code embodied in the medium, the computer 
readable program code configured to carry out the method of 
any of claims 14 to 27 or 66 to 92 of the appended claims. 
[0045] An advantage of the invention is that effective 
detection of threat objects at stand-off distances maybe 
accomplished. 

[0046] A further advantage of the invention is that detec- 
tion of threat objects using a portable device may be 
achieved. 

[0047] A further advantage is that the detection of smaller 
dimensions (e.g. gun barrel dimensions and/or caliber) is 
enabled. 


BRIEF DESCRIPTION OF THE FIGURES 


[0048] Embodiments of the invention will now be 
described in detail, by way of example, with reference to the 
accompanying drawings, in which: 

[0049] FIG. 1 is a block diagram of an object detection 
system according to some embodiments of the invention; 
[0050] FIG. 2A shows a block diagram of an object 
detection system according to some embodiments of the 
invention. 

[0051] FIGS. 2B and 2C show the different behaviour of 
smooth and jagged conducting targets in response to linearly 
polarized microwave radiation; 

[0052] FIG. 3 shows the original data (amplitude vs fre- 
quency) for the radiation detected by the system of FIG. 1 
from the scan of an x-band waveguide of major dimension 
80 mm; 

[0053] FIGS. 4A and 4B illustrate, respectively, Burg and 
Fourier transforms of the data set of FIG. 3; 

[0054] FIG. 5 shows a sample metal object (an X band 
waveguide connector) with dimensions, as used in testing 
the system of FIG. 1 to generate the data of FIG. 3; 
[0055] FIG. 6 shows a comparative plot of the Fourier 
Transformed data taken from a repeated rapid scan between 
14 to 40 GHz at a standoff distance of a person with small 
handgun that was strapped to the front of the person, as the 
person rotated round slowly within the beam, using the 
system of FIG. 1 or 2A, and using swept reflectrometry 
techniques according to embodiments of the invention; 
[0056] FIG. 7 illustrates the same plot as in FIG. 6, but 
with the handgun absent; 

[0057] | FIG. 8A shows a plot using the Burg transform for 
a scan of a person with simulated plastic explosive strapped 
to their midriff, using swept reflectrometry techniques 
according to embodiments of the invention. 
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[0058] FIG. 8B shows the same plot as in FIG. 8A, but 
with the simulated plastic explosive absent; 

[0059] FIG. 9 illustrates measured frequency response of 
detected signal from a metal barrel with axis aligned at the 
angles shown to the microwave propagation direction, using 
barrel tone detection techniques according to some embodi- 
ments of the invention; 

[0060] FIG. 10 shows correlation of the observed signal 
with the calculated chirped response function for three 
different values of fo and the direct Fourier transform, using 
barrel tone detection techniques according to some embodi- 
ments of the invention; 

[0061] FIG. 11 shows representative signals, displayed in 
the time domain, from a number of targets—in each case the 
signal has been measured simultaneously with separate 
horns in two polarizations, using cross-polarization detec- 
tion techniques according to some embodiments of the 
invention; 

[0062] FIG. 12 illustrates the scanning of a simple metal 
plate using later time response techniques according to some 
embodiments of the invention, with the EM field polarized 
in the direction of the 5 cm dimension; 

[0063] FIG. 13 shows the time domain response of the 5 
cm wide flat sheet of FIG. 12, rotated through various angles 
between 0 and 60 degrees, using later time response tech- 
niques according to some embodiments of the invention; 
[0064] FIG. 14 shows a typical late time responses from a 
concealed handgun, including the cross-polarized signals 
(crossed line), using later time response techniques accord- 
ing to embodiments of the invention; 

[0065] FIG. 15 shows schematically the use of a neural 
network to make decisions on whether a threat object is 
present, based on inputs determined by techniques according 
to some embodiments of the invention; 

[0066] FIG. 16A shows the arrangement for detection of 
the Late Time Response according to a further embodiment 
of the invention in plan view. 

[0067] FIG. 16B shows the arrangement for detection of 
the Late Time Response as viewed from A-A cross-section. 
[0068] FIG. 17 shows a plan view for a detection system 
according to a further embodiment of the invention incor- 
porating multiple transmitter and receiver pairs, 

[0069] FIGS. 18A and 18C show a typical response of a 
human body with concealed handgun, FIG. 18B shows a 
similar trace with the handgun absent; 

[0070] FIGS. 19A, 19B, and 19C show the process of 
detecting objects using LTRs in accordance with embodi- 
ments of the invention; 

[0071] FIGS. 20A and 20B show the typical LTR response 
for representative metal objects 

[0072] FIG. 20C shows the typical LTR response for a 
complex object (FIG. 20C) such as a handgun; 

[0073] FIG. 21 illustrates the total time response for a 
Glock 17 handgun suspended in air, showing the LTR (3), to 
be separated from the other components of the signal; and 
[0074] FIGS. 22A, 22B and 22C illustrate poles extracted 
from LTR data using General Pencil of Functions method, 
for a spanner, an Allen key, and a handgun, respectively. 


DETAILED DESCRIPTION OF EMBODIMENTS 


[0075] Embodiments of the present invention now will be 
described more fully hereinafter with reference to the 
accompanying drawings, in which embodiments of the 
invention are shown. This invention may, however, be 
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embodied in many different forms and should not be con- 
strued as limited to the embodiments set forth herein. Rather, 
these embodiments are provided so that this disclosure will 
be thorough and complete, and will fully convey the scope 
of the invention to those skilled in the art. Like numbers 
refer to like elements throughout. 

[0076] It will be understood that, although the terms first, 
second, etc. maybe used herein to describe various elements, 
these elements should not be limited by these terms. These 
terms are only used to distinguish one element from another. 
For example, a first element could be termed a second 
element, and, similarly, a second element could be termed a 
first element, without departing from the scope of the present 
invention, As used herein, the term “and/or” includes any 
and all combinations of one or more of the associated listed 
items. 

[0077] It will be understood that when an element such as 
a layer, region or substrate is referred to as being “on” or 
extending “onto” another element, it can be directly on or 
extend directly onto the other element or intervening ele- 
ments may also be present. In contrast, when an element is 
referred to as being “directly on” or extending “directly 
onto” another element, there are no intervening elements 
present. It will also be understood that when an element is 
referred to as being “connected” or “coupled” to another 
element, it can be directly connected or coupled to the other 
element or intervening elements may be present. In contrast, 
when an element is referred to as being “directly connected” 
or “directly coupled” to another element, there are no 
intervening elements present. 

[0078] Relative terms such as “below” or “above” or 
“upper” or “lower” or “horizontal” or “vertical” may be used 
herein to describe a relationship of one element, layer or 
region to another element, layer or region as illustrated in the 
figures. It will be understood that these terms are intended to 
encompass different orientations of the device in addition to 
the orientation depicted in the figures. 

[0079] The terminology used herein is for the purpose of 
describing particular embodiments only and is not intended 
to be limiting of the invention. As used herein, the singular 
forms, “a”, “an” and “the” are intended to include the plural 
forms as well, unless the context clearly indicates otherwise. 
It will be further understood that the terms “comprises” 
“comprising,” “includes” and/or “including” when used 
herein, specify the presence of stated features, integers, 
steps, operations, elements, and/or components, but do not 
preclude the presence or addition of one or more other 
features, integers, steps, operations, elements, components, 
and/or groups thereof. 

[0080] Unless otherwise defined, all terms (including tech- 
nical and scientific terms) used herein have the same mean- 
ing as commonly understood by one of ordinary skill in the 
art to which this invention belongs. It will be further 
understood that terms used herein should be interpreted as 
having a meaning that is consistent with their meaning in the 
context of this specification and the relevant art and will not 
be interpreted in an idealized or overly formal sense unless 
expressly so defined herein. 

[0081] As used herein, “threat object” is taken to mean a 
metallic or dielectric object, whether specifically designed 
or intended for offensive use or not, that have potential to be 
used in an offensive or violent manner. 

[0082] The present invention is described below with 
reference to flowchart illustrations and/or block diagrams of 
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methods, systems and computer program products according 
to embodiments of the invention. It will be understood that 
some blocks of the flowchart illustrations and/or block 
diagrams, and combinations of some blocks in the flowchart 
illustrations and/or block diagrams, can be implemented by 
computer program instructions. These computer program 
instructions may be stored or implemented in a microcon- 
troller, microprocessor, digital signal processor (DSP), field 
programmable gate array (FPGA), a state machine, pro- 
grammable logic controller (PLC) or other processing cir- 
cuit, general purpose computer, special purpose computer, or 
other programmable data processing apparatus such as to 
produce a machine, such that the instructions, which execute 
via the processor of the computer or other programmable 
data processing apparatus, create means for implementing 
the functions/acts specified in the flowchart and/or block 
diagram block or blocks. 

[0083] These computer program instructions may also be 
stored in a computer readable memory that can direct a 
computer or other programmable data processing apparatus 
to function in a particular manner, such that the instructions 
stored in the computer readable memory produce an article 
of manufacture including instruction means which imple- 
ment the function/act specified in the flowchart and/or block 
diagram block or blocks. 

[0084] The computer program instructions may also be 
loaded onto a computer or other programmable data pro- 
cessing apparatus to cause a series of operational steps to be 
performed on the computer or other programmable appara- 
tus to produce a computer implemented process such that the 
instructions which execute on the computer or other pro- 
grammable apparatus provide steps for implementing the 
functions/acts specified in the flowchart and/or block dia- 
gram block or blocks. It is to be understood that the 
functions/acts noted in the blocks may occur out of the order 
noted in the operational illustrations. For example, two 
blocks shown in succession may in fact be executed sub- 
stantially concurrently or the blocks may sometimes be 
executed in the reverse order, depending upon the function- 
ality/acts involved. Although some of the diagrams include 
arrows on communication paths to show a primary direction 
of communication, it is to be understood that communica- 
tion may occur in the opposite direction to the depicted 
arrows, 

[0085] Embodiments of the invention may be used for 
remotely detecting the presence and/or size of metal and/or 
dielectric objects concealed underneath clothing. Embodi- 
ments herein may be used for remotely detecting metal 
and/or dielectric objects. A dielectric in this context is a 
non-conducting (i.e. insulating) substance such as ceramic 
that has a low enough permittivity to allow microwaves to 
pass through. A ceramic knife or gun, or a block of plastic 
explosive, are examples of this type of material. 

[0086] Some embodiments of detection systems are dis- 
closed herein, FIG. 1 includes embodiments using direct 
detection without phase detection and FIG. 2 includes phase 
detection embodiments. In some embodiments, the hard- 
ware may be embodied in a portable and covertly deployable 
system. 

[0087] FIG. 1 is a block diagram of the threat object 
detection system according to some embodiments of the 
invention. For the direct detection (without the phase) 
responses, a detection system may include a microwave 
and/or mm wave source 102 (40 GHz Agilent Microwave 
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Synthesizer). In some embodiments, the microwave and/or 
mm wave source 102 may comprise a separate component of 
the system. In some embodiments, the detection system may 
include a controller (PC) 104, two 20 dB standard gain horns 
used as transmitter 106 and receiver 108 for the K.0 and Q 
bands, a zero-bias direct detector 110 followed by a DC 
amplifier 112, and a high speed data acquisition card (PCI- 
6132 National Instalment interface). In some embodiments, 
the system maybe controlled using control software includ- 
ing Labview or C# code, among others. In some embodi- 
ments, free wave directional antennae may also replace 
horns to permit a wider scanning range. 


[0088] The main procedure carried out by controller 104 
in implementing a scan is set out in Appendix A at the end 
of this section of the Specification. The procedure “Perform 
transformation on received radiation signals to produce time 
domain or optical depth domain trace” in Appendix A, for 
each of the four techniques, is discussed later with reference 
to Appendix B. 


[0089] In use and operation, the system may use electro- 
magnetic radiation in the microwave or millimetre (mm) 
wave band, where the wavelength is comparable or shorter 
than the size of the object 116 to be detected. The object 116 
may be on and/or in the body of a person, within containers 
and/or items of luggage, and/or concealed in and/or on some 
other entity (not shown). The suspect entity (e.g., a person; 
not shown) has radiation directed by transmitter 106 onto it, 
so that the (threat) object 116 is entirely illuminated by a 
continuous wave of this radiation (i.e., the radiation is not 
pulsed, but kept continuously on). The radiation intensity is 
well within safe operating limits, but may be in any case 
determined by the sensitivity of the detector 110, As an 
example, in the range 14-40 GHz, 0 dBm of power is used 
with a typical beam area 118 of 0.125 m?, which equates to 
a 20 cm diameter beam. However, in some embodiments, the 
hardware maybe designed so as generate a beam area 118 of 
greater or lesser size. 


[0090] The frequency and consequently the wavelength of 
the radiation, is swept through a reasonable range and may 
be referred to as swept CW and/or continuous wave radia- 
tion. Limits maybe set by the devices used, but a 20 GHz or 
more sweep starting at 14, 50 or 75 GHz is typical. The data 
may be collected in a series of small steps and/or as a 
real-time continuous sweep. Typically 256 or more data 
points may be acquired. In some embodiments, data may be 
taken between 14 to 40 GHz, providing a sweep range of 26 
GHz. 


[0091] The illumination and detection may be undertaken 
remotely from the object 116 in question, for example, at a 
distance of a meter or more, although there is no lower or 
upper limit on this distance. The upper limit on detection 
distance may be set by the millimeter or microwave focusing 
optics, although, with this technique, a small beam at the 
diffraction limit is not necessary. In some embodiments, 
ranges for this device may include a few tens of centimetres 
(cm) to many tens of meters (m), In some embodiments, a 
device may be operated at a range of approximately 1 m to 
10 m depending on the frequency chosen (some microwave 
frequencies are attenuated by the atmosphere, and atmo- 
spheric windows such as that found around 94 GHz are 
generally chosen to minimise these effects) In some 
embodiments, the source of electromagnetic radiation 102 
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and the detector 110 may be mounted next to each other and 
they may be focused onto some distant object 116 or entity 
(not shown). 


[0092] A variety of techniques is disclosed herein and may 
include distinct system and/or method recitations. For 
example, swept reflectrometry and barrel tone detection, 
may provide that the return radiation is detected and its 
amplitude stored as a function of frequency. In this regard, 
embodiments of a system as illustrated in FIG. 1 may be 
used. 


[0093] Other techniques, including cross-polarization and 
LTR recurrences may use the phase of the returned radiation 
that may be acquired at each frequency point. When the 
phase of the returned signal is used—in order to replicate a 
Time Domain response via the use of a Fourier Transform— 
the synthesiser and detection system in FIG. 1 may be 
replaced by a four port Vector Network Analyser (VNA) (for 
example 40 GHz Rohde Schwarz SVA Vector Network 
Analyzers) or other device capable of obtaining the phase 
and amplitude of the returned signal, once again swept from 
14 to 40 GHz in small steps. 


[0094] FIG. 2A shows a block diagram of an object 
detection system according to some embodiments of the 
invention. As illustrated, components of the VNA may be 
included in some embodiments, Others of the system com- 
ponents are discussed above regarding FIG. 1, except as 
described below. 


[0095] The frequencies of the first and second microwave 
sources 102 and 103 may be swept under control of the 
signal from the Ramp Generator 202 to remain approxi- 
mately 100 MHz apart. The Microwave mixers 204, 206 
generate signals corresponding to the difference frequency 
between the two inputs (-100 MHz). After amplification by 
two RF amplifiers 208, 210, a RF Mixer 212 produces two 
outputs corresponding to the “in phase" (T) and “in quadra- 
ture" (Q) components of the detected signal. The signals 
may be amplified by amplifiers 112, 112', and the data 
acquisition may be controlled by a controller 104 (PC). The 
entire system apart from the horns 106, 108 may be referred 
to as a Vector Network Analyser (VNA). 


[0096] The return signal is collected by a horn and applied 
to port 2 of the VNA, which measures parameter S21, If 
cross-polarization measurements are used, a second receiver 
horn oriented (not shown) at 90° may be added to the VNA 
on, for example, port 3, The transmitted signal may be 
generated from port 1 and may be, for example, 1 mW. The 
real and imaginary parts may be recorded and can be 
corrected for the electrical behaviour of the horns. The 
signals are zero padded out to 4096 points and are processed 
by a Fast Fourier Transform routine to yield the effective 
time response. In this manner, the process allows the repli- 
cation of the application of a pulse of radiation to the target 
(entity) and the subsequent acquisition of the time resolved 
response. 


[0097] FIGS. 2B and 2C show the different behaviour of 
smooth and jagged conducting targets in response to linearly 
polarized microwave radiation. A single reflection from a 
metal surface produces a signal polarized in the same 
manner as the incident signal—polarization is conserved as 
illustrated in FIG. 2B, However, multiple reflections 
on/around the target will rotate the direction of polarization 
as illustrated in FIG. 2C. In some embodiments, this can be 
detected by using an additional waveguide horn (not shown) 
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that is rotated 90 degrees relative to the first horn 108, thus 
blocking the polarization conserving signal. 

[0098] In some embodiments, hardware corresponding to 
the systems herein may form and/or be part of a portable 
device (i.e. small enough to be carried by one person, or 
transported in an automobile, so as to be operable therein), 


Theoretical Basis. 


[0099] To detect the range of signals described herein, 
several measurement techniques are available, For a trans- 
mitted signal the return signal Ex. from a target (entity) 
distance L away may be written as follows: 


Ep=rEoe te? “f° (1) 


where co-2nf, f is the frequency, r the scattering coefficient 
and c the velocity of light. 

[0100] A detector 110, which may respond to the micro- 
wave power, will only measure |E,|?, which for a single 
scattering center, does not explicitly depend on frequency. 
However, for two scattering segments at different ranges L,, 
and L,, the power is proportional to: 


2344. 242 |? (2) 
IER? = nie € one c | 


[0101] This contains terms in cos (2w(L,)c) (i.e., oscilla- 
tory terms) that are dependent on the difference in range 
L,-L,. By performing a Fourier transform on the detected 
power measured as a function of transmitted frequency, 
peaks corresponding to the difference in range of various 
parts of the target are observed and these give an indication 
of the size of the object 116. However, for a complicated 
object 116, many pairs of distances would be involved and 
the analysis of the signal would be complex. 

[0102] For a different group of detectors (i.e. Vector Net- 
work Analysers as illustrated in FIG. 2A), it is possible to 
measure the complex return signal directly, These Vector 
Network Analysers effectively mix the return signal with a 
fraction of the transmitted wave to measure 


2.0] 21.92 
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in terms of its real and imaginary parts. In this case Fourier 
Transforming leads to a series of signal peaks at the range 
of each element of the target and arranged in the order of 
their distance. Thus, a much clearer indication of the dimen- 
sions of the object 116 may be obtained, though only in one 
dimension. Any Late Time Responses (i.e. those that cannot 
be attributed to direct scattering) can be measured in this 
way, although their strength may be many times less than 
directly reflected signals. 

[0103] Further information about the target (entity) may 
be obtained using a second detector to collect return signals 
emitted at a different angle from the target. This effectively 
probes the target along a second direction and can in 
principle enable more dimensions of the object 116 to be 
ascertained. 

[0104] A Fourier Transform (FT) or some other more 
advanced power spectrum analysis technique such as a Burg 
Transform may be applied. The Burg and related methods of 
power spectrum analysis may be better than the FT for this 
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application as the individual peaks that relate directly to the 
dimensions of the object are more clearly identifiable, and as 
it is possible to choose the number of peaks to be displayed 
in the output (and hence reject weaker peaks), They may also 
allow two closely spaced peaks to be resolved. 


[0105] FIG. 3 shows original experimental data (ampli- 
tude vs frequency) for the radiation detected by the system 
of FIG. 1 from the scan of an x-band waveguide of major 
dimension 80 mm. FIGS. 4A and 4B illustrate respectively 
Burg and Fourier transforms of the data set of FIG. 3. FIG. 
5 shows a sample metal object (an X band waveguide 
connector) with dimensions, as used in testing the system of 
FIG. 1 to generate the data of FIG. 3. The Fourier Transform 
(FIG. 4B) and Burg Transform (FIG. 4A) are presented for 
comparison. The peak at 80 mm corresponds with the length 
of the waveguide. The Burg algorithm as identified in FIG. 
4A is much less cluttered than a conventional FT as illus- 
trated in FIG. 4C. The Burg algorithm is used here to turn 
the frequency sweep into a power spectrum, In this regard, 
the peaks in the power spectrum may directly give the 
various dimensions and/or lengths of the metal and/or 
dielectric object, after appropriate scaling. 

[0106] The position L of peaks within the FFT or Burg 
spectrum directly relate to the size of the object using the 
following formula: 


c 
287 


Where c is the speed of light, Af the periodicity in the 
freguency domain. This azis represents optical depth for the 
purposes of this disclosure. 

[0107] The minimum spatial resolution d is related to the 
sweep range or bandwidth BW: 


2BW 


[0108] As an example, if the source frequency were to be 
swept between 14 and 40 GHz, this constitutes a sweep 
range of 26 GHz, which translates to a resolution of 5.7 mm. 
A larger sweep range would lead to an improved resolution, 
which may result in, for example, a maximum optical depth 
or distance of 740 mm for 256 data points. In some embodi- 
ments, the number of data points may include more than 256 
including, for example, 512, 1024 or any multiple thereof, 
among others, For Vector Network Analysers operating in 
Time Domain mode, in which the complex data is converted 
to the time domain, this calculation may be included in the 
software, 

[0109] The four techniques mentioned above, each of 


which may be used in embodiments of the invention, will 
now be discussed in more detail. 


Technique 1: Swept Reflectometry. 


[0110] As briefly described above, swept reflectometry is 
the principle by which the distances between corners, edges 
and cavities on the threat object (weapon) 116 independently 
of the distances to the source (TX horn 106; FIG. 1) and the 
detector (RX horn 108; FIG. 1) maybe ascertained in real- 
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time. Here, 11312 as described above is the measured quan- 
tity. The distance to the target (entity) may not always be 
available. 


[0111] If the object 116 (e.g. a weapon strapped to the 
body as the latter rotates) is moved around in the beam and 
its angle and distance with respect to the source and detector 
is changed, those dimensions between edges and corners 
that actually belong to the object can be differentiated from 
those that do not. In this manner, the background clutter may 
be removed. 


[0112] Some embodiments of a procedure carried out by 
controller 104 (FIG. 1) in implementing a scan is set out in 
Appendix A at the end of this section of the disclosure. The 
procedure “Perform transformation on received radiation 
signals to produce time domain or optical depth domain 
trace" for Technique 1 is set out in Appendix B.1 at the end 
of this section of the disclosure. 


[0113] The software according to Technique 1 may differ- 
entiate the peaks that relate to the dimensions of the object 
116 from those that do not by acquisition of the signal over 
a period of time and by storing these acquired signals 
independently, with the object moving within the beam. The 
signals that indirectly relate to the dimensions of the object 
remain and occur within certain bands denoted by the 
distances between the various corners of the object. Other 
signals that change (e.g., the air gaps between clothing and 
the skin) may be more chaotic and may be integrated out 
over a period of time. 


[0114] Ifthe strength ofthe signal is normalized relative to 
distance, the returns from a subject concealing a handgun 
will be larger in amplitude than those without. 


[0115] FIG. 6 shows a comparative plot of the Fourier 
Transformed data taken from a repeated rapid scan between 
14 to 40 GHz at a standoff distance of a person with small 
handgun strapped to the front of the person, as the person 
rotated round slowly within the beam, using swept reflec- 
trometry techniques according to embodiments described 
herein. FIG. 7 illustrates the same plot as in FIG. 6, but with 
the handgun absent. The scans are presented in three dimen- 
sions, with the scan number on the X axis, optical depth on 
the Y axis and the amplitude of the power spectrum (arbi- 
trary-units) on the Z axis. In FIG. 6, many resonances can be 
seen below 100 mm with the gun present, denoting the 
various distances between corners and edges of the weapon. 
The second plot (FIG. 7) shows the response from the body 
alone at the same standoff distance and with the body 
rotating in the same manner. 


[0116] Very large dimensions, such as, for example, metal 
doors, window frames and/or a multitude of other metal 
objects will not be observed as they are not entirely encom- 
passed by the beam, as the microwave beam can be focused 
onto relevant parts of the person (entity) or object 116 in 
question. 


[0117] The reflected return radiation is seen to contain 
patterns or frequencies that can be indirectly related to the 
dimensions of the metal object according to the technique 
identified in the section “Theoretical basis" above, including 
the presence and/or length of gun barrels. In this manner, 
embodiments herein may be used to discriminate between, 
say, handguns and keys, knives and keys, etc. In effect, the 
technique measures the distances between the various edges 
ofthe object at the orientation of the source and detector, and 
cavity lengths if present. 
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[0118] The dimensions of guns and knives are different 
from most other objects carried about the person, so the 
appropriate dimensions may be stored on a database. 
[0119] The technique (Technique 1) is also capable of 
measuring, particularly, the depth of a dielectric (i.e., of a 
material that does not conduct electricity), although the 
physics behind this may be significantly different, A dielec- 
tric object might typically be a lump of plastic explosives 
concealed on a suicide bomber. 

[0120] FIG. 8A shows a plot using the Burg transform for 
a scan of a person with simulated plastic explosive strapped 
to their midriff, using swept reflectrometry techniques 
according to some embodiments of the invention, and FIG. 
8B shows the same plot as in FIG. 8A, but with the simulated 
plastic explosive absent. The Burg transform (FIG. 8A) is of 
a 14-40 GHz scan of a person carrying an 80 mm thick block 
of plastic explosive simulant with a dielectric constant of 
1.5, giving an apparent optical depth of 120 mm. 


Technique 2: Barrel Tone Detection by Direct Detection of 
Aspect-Independent Chirped Signals. 


[0121] Threat objects that contain cavities and can be 
excited by an incoming microwave signal, can exhibit strong 
frequency dependence in the scattered signal. These signals 
may differ from those derived from the outside of the object 
by: 


[0122] 1. Showing a threshold frequency for stimulation 
(cut-off), 
[0123] 2. Being less dependent on alignment of the 


cavity with respect to the microwave direction. 


[0124] For example, consider a 10 cm long cylindrical 
metal barrel closed at one end, which has 19 mm outside 
diameter and 9 mm inside diameter, The HI 1 mode has the 
lowest threshold frequency/o for propagation for the inside 
bore at 19.5 GHz. For a cavity length L, the response 
for/greater than/o is proportional to the chirped sine wave 
signal: 


IEP ৭৩০21577102] ৰ 


where/is the microwave freguency and c is the velocity of 
light, with a minimal return at lower freguencies below this 
threshold. 


[0125] A procedure according to some embodiments that 
may be carried out by controller 104 (FIG. 1) in implement- 
ing a scan is set out in Appendix A at the end of this section 
of the disclosure, The procedure “Perform transformation on 
received radiation signals to produce time domain or optical 
depth domain trace” for Technique 2 is set out in Appendix 
B.2 at the end of this section of the disclosure. 


[0126] FIG. 9 illustrates measured frequency response of 
detected signal from a metal barrel with axis aligned at the 
angles shown to the microwave propagation direction, using 
barrel tone detection techniques according to the embodi- 
ments of the invention. As shown in FIG. 9, the return signal 
clearly shows the onset of the oscillatory response above the 
threshold/o. In this regard, the onset of chirped oscillations 
occurs when the microwave frequency is beyond cut-off for 
propagation through the bore of the barrel. The bottom trace 
shows the optimum calculated chirped response. 
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[0127] It should be noted that the signal is clearly seen at 
a range of impact angles @ ranging from 0 to over 45° and 
the oscillation frequency is only a function of the cavity 
length L. This may be contrasted with the case of edge or 
corner detection where the oscillation frequency is propor- 
tional to L cos(f). The analysis provides a means of deter- 
mining both the length and diameter of the cavity bore. 
[0128] FIG. 10 shows correlation of the observed signal 
with the calculated chirped response function cos(2n(f?—fo*) 
“2(2L/c)+@) for three different values of fo and the direct 
Fourier transform (fo=0) using barrel tone detection tech- 
niques according to some embodiments of the invention. In 
FIG. 10 the measured signal is correlated with the chirped 
wave response for particular values of threshold frequency 
fo. The correlated signal clearly peaks at the actual length L, 
i.e. 10 cm. It can be seen that this signal is much sharper than 
the direct FT, i.e. when fo=0 and the signal is strongly 
dispersed. Stated differently, the correlation function is 
sharper and more symmetric when/o is equal to the true 
threshold frequency of 19500 MHz. It is also significantly 
narrower and more symmetric than when the assumed value 
for fo is higher or lower than the correct value. This provides 
a method of determining both the diameter and length of a 
cylindrical cavity and thus identifying a threat object such as 
a gun barrel. 


Technique 3: Identification of the Target by Cross-Polarized 
Detectors. 


[0129] By measuring the return signal as a function of 
frequency scanned over a wide range, the dimensions may 
be recovered through Fourier Transform techniques. This 
duplicates the effect of responses of targets, to a very short 
excitation pulse without the need of high speed switches and 
ultra-fast detectors and digitization processes. The range 
resolution obtainable is of the order of 0.5-1.0 cm, as 
described by the principles above, at the sweep ranges 
available here (14-40 GHz, but this property is not restricted 
to this frequency range), and thus appropriate for character- 
izing objects such as hand guns. 

[0130] Another useful aid to threat object identification is 
to measure the polarization of the return signal. Waveguide 
horns (see FIG. 2A) act as excellent polarizers and if the 
transmit 106 and receive 108 horns are similarly oriented, 
then polarization conserving components are detected. How- 
ever, a second horn rotated about its axis by 90 degrees will 
be blind to “normal” polarization conserving signals and 
only detect “cross” or polarization changing signals. Con- 
ducting materials with a smooth surface, including the 
human body, are mainly polarization conserving, However, 
complicated targets that involve multiple reflections at dif- 
ferent angles and particularly metal objects with sharp 
edges, give rise to significant “cross” polarized signals, 
which can lead to good discrimination, for example of a 
hand gun next to the body. 

[0131] Embodiments of a procedure carried out by con- 
troller 104 (FIG. 2A) in implementing a scan is set out in 
Appendix A at the end of this section of the disclosure. The 
procedure “Perform transformation on received radiation 
signals to produce time domain or optical depth domain 
trace” for Technique 3 is set out in Appendix B.3 at the end 
of this section of the disclosure. 

[0132] The responses of a range of objects using a system 
as illustrated in FIG. 2A with a second receiver horn oriented 
(not shown) at 90° and with a cross-polarized detector (not 
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shown) placed almost immediately above the first detector 
are analyzed. The test targets ranged from a human alone, 
front and side, the same configuration with a concealed 
small handgun, a bunch of keys, a mobile phone and a digital 
camera. 


[0133] FIG. 11 shows representative signals, displayed in 
the time domain, from a number of targets—in each case the 
signal has been measured simultaneously with separate 
horns in two polarizations, using cross-polarization detec- 
tion techniques according to embodiments of the invention. 
FIG. 11 shows the relevant section of the responses, illus- 
trating distance and object information before and after the 
target has been removed. In each case, the signal has been 
measured simultaneously with separate horns in two polar- 
izations, parallel to the transmitted beam (single line) and at 
right angles (crossed line). The quoted range is only relative, 
in some embodiments the targets are 1-2m away from the 
horns, Signals (a) and (b) are from a small hand gun, (c) and 
(d) from the chest area of a human body in different 
orientations, (e) and (f) are for the gun held next to the body, 
(g) for a small camera held, and (h) for a set of keys. It can 
be can be seen that quite distinctive behaviour is found for 
the gun in the two polarizations when compared to the body 
alone and with objects with flat surfaces. 


[0134] In accordance with Technique 3, the very wide 
sweep range leads to detailed information about the dimen- 
sions of the object 116 and the fine structural distances 
between the source/detector 106/108 and corners/edges on 
the target/object 116. It can be seen in FIG. 11A that the 
signal from the cross-polarized detector (crossed, line) is not 
only enhanced by the presence of the weapon, but the 
structural detail is to some extent anti-correlated. Compare 
this plot to the one taken from the body-only, where the cross 
polarized detector shows little information. Similarly, when 
the gun is placed on the body (either in front or at the side), 
the cross-polarized signal is once again enhanced, with a 
degree of anti-correlation between the normal (solid line) 
and cross polarized (crossed line) detectors. 


Technique 4: Target Determination by Aspect Independent 
Effects (Late Time Responses). 


[0135] The Late Time Response (LTR) and the closely 
associated Singularity Expansion method (SEM) arose from 
the observation that the time-resolved radar signature from 
conducting objects contains information after the radar pulse 
has passed the target. 


[0136] The pulse sets up currents on the surface of the 
object 116 in the form of resonant modes, which subse- 
quently re-radiate. An alternative interpretation is to con- 
sider the radar pulse stimulating travelling waves on the 
surface of the target, which move across and around the 
object 116 until they return back to their initial distribution, 
This recurrence can re-radiate back into the return beam, 
which appears an extra, time-delayed signal. 


[0137] An important feature of the LTR is that the time 
taken for the travelling wave to circulate around the object 
does not depend on the orientation of the entity/object 116, 
and hence is aspect-independent, although its strength may 
depend on the efficiency of coupling into the modes. For 
objects 116 with symmetry the time-delay is approximately 
half the perimeter with respect to centre of the object 116. 
For more complicated objects, the LTR signal is more 
complex, but the structure can be interpreted in terms of the 
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dimensions of the object 116, These include dimensions 
cross range as well as the usual along range values. 
[0138] FIG. 12 illustrates the scanning of a simple metal 
plate 1202(PCB) using later time response techniques 
according to some embodiments of the invention. The return 
signal from a rectangular piece of copper coated PCB was 
measured. Its narrow dimension (50 mm) is oriented in the 
plane of the microwave electric field and the long dimension 
(200 mm) is at right angles to the polarization direction and 
the direction of propagation. 

[0139] The return signal is collected by a second horn 
close to the transmitting horn and is measured and applied, 
for example, to port 2 of a VNA 1204. The S21 parameter 
is recorded over a range typical range 14-40 GHz. The 
signals were corrected for the performance ofthe microwave 
horns 106, 108. 

[0140] If the plate 1202 of width L is rotated by angle 6 
about its long direction, then scattering from its leading and 
trailing edge leads to a doublet response with separation L 
sin 0. 

[0141] A procedure according to some embodiments may 
be carried out by controller 104 (FIG. 2A) in implementing 
a scan and is set out in Appendix A at the end of this section 
of the disclosure. The procedure “Perform transformation on 
received radiation signals to produce time domain or optical 
depth domain trace” for Technique 4 is set out in Appendix 
B.4 at the end of this section of the disclosure. 

[0142] FIG. 13 shows the time domain response of the 5 
cm wide flat sheet of FIG. 12, rotated through various angles 
between 0 and 60 degrees, using later time response tech- 
niques according to some embodiments of the invention. 
The scattering from leading and trailing edges is centered 
around 10 cm. There is also a clear Late Time Response at 
15 cm whose amplitude is largely independent of angle. Also 
seen in FIG. 13 is a LTR response at a fixed distance of 5 cm. 
[0143] FIG. 14 shows late time responses from a con- 
cealed handgun, including the cross-polarized signals 
(crossed line), using later time response techniques accord- 
ing to some embodiments of the invention. This shows the 
decaying oscillatory responses after 30 cm, being typical of 
late time responses from complex objects. The cross-polar- 
ized signal is denoted by the crossed line. 

[0144] FIG. 15 shows schematically the use of a neural 
network 1500 to make decisions on whether a tin-eat object 
is present, based on inputs 1502 determined by techniques 
according to some embodiments of the invention. The 
outputs 1504 from the one or more of the abovementioned 
Techniques 1 to 4 may be taken to the neural network 
classifier (e.g., aback propagation feed forward network) 
that has been pre-trained on sets of the type of concealed 
threat objects 116 that will be of concern, in addition to 
harmless items. 

[0145] The training data may be formed of sets of data 
taken using the above-described techniques, in random 
order, The output 1506 from the neural network may be a 
single output that gives a confidence level (1=gun, 0=no 
object being concealed). However, it will be appreciated that 
other configurations or learning algorithms may be used. For 
example, multiple outputs from the neural network 1500 
maybe employed, with sub-classifications (e.g., gun, mobile 
phone, keys, etc). 

[0146] Tine millimeter wave reflected signals from the 
guns show a number of features which enable them to be 
distinguished from innocently carried objects such as keys 
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and mobile phones. These include cavity mode oscillations 
from the barrel with characteristic frequency of onset that 
allows the caliber and length to be “determined. The inter- 
ference of signals from different parts of the target leads to 
a frequency-dependent response, which can be used to 
deduce the size of the object. The response at different 
polarizations may give an indication of the complexity of the 
object from multiple reflections. The late time response 
gives an aspect-independent signal dependent on target 
dimensions. Taken together these features provide a means 
of detecting handguns, for example, under practical condi- 
tions at stand-off distances. The application of signal pro- 
cessing techniques enable relevant parameters to be 
extracted for use in automatic detection systems. 

[0147] Embodiments of the system forming extensions 
and variants of Technique 4: Target Determination by aspect 
independent effects (Late Time Responses), will now be 
discussed in more detail. 

[0148] FIG. 16 shows the arrangement for detection using 
LTRs, according to a further embodiment of the invention, 
(a) in plan view, and (b) as viewed at A-A in FIG. 16A. A 
TX horn 106' and a RX horn 108' are mounted facing the 
space in which an entity or target (e.g. person) 116, here 
shown carrying a metal threat object 1602, stands or may 
pass, during scanning for threat objects. Preferably, TX horn 
106', and TX horn 108', each comprise a dual polarised horn 
antenna. 

[0149] The TX horn 106' is coupled to detection electron- 
ics 1604 via (coax) cables 1606, and RX horn 108' is 
coupled to detection electronics 1604 via (coax) cables 
1608. The detection electronics 1604 may comprise the 
circuitry of FIG. 2A, excluding the horns 106,108 and may, 
to the extent required, including PCs, processing circuitry 
and software controlled devices, for implementing the scan- 
ning, detection and data analysis techniques according to 
this embodiment, as discussed earlier in relation to other 
embodiments. 


[0150] In certain embodiments, the apparatus includes a 
linear polarising element 1610, mounted for rotation about 
axis B. In this way, the entity or target 116 is illuminated by 
electromagnetic field with rotating polarisation. The spin- 
ning linear polarising element 1610 is used in conjunction 
with a dual polarised horn antenna to produce broadband 
microwaves with polarisation which rotates (either ellipti- 
cally or circularly, depending on the relative amplitudes of 
the dual microwave feed) at a frequency equal to that of the 
frequency of rotation of the linear polarising element 1610. 
This allows optimum coupling between the exciting wave 
and the target 116, since coupling is strongly dependent 
upon polarisation for a given target aspect. The receiver RX 
horn 108' is also a dual polarised antenna and data from both 
channels are used in processing, as discussed further here- 
inafter. 

[0151] FIG. 17 is a plan view of a configuration for a 
detection system according to a further embodiment of the 
invention, incorporating multiple transmitter and receiver 
pairs 1702, Suitably, each transmitter and receiver pair 1702 
comprises the TX horn 106' and a RX horn 108' of FIG. 16. 
In this case, four transmitter and receiver pairs 1702 are 
used, positioned in an anechoic lined corridor 1704 to give 
all round interrogation as the person walks through (in the 
path approximately indicated by arrows C). The use of 
multiple transmitter and receiving antennae pairs 1702 func- 
tions so as to give optimal coverage as the person (not 
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shown) being interrogated passes through a corridor 704. 
These antennae pairs 1702 are concealed behind microwave 
transparent windows 1706 in walls 1708 rendered approxi- 
mately anechoic by use of microwave absorbing foams 
1710. The microwave receiver antenna is a dual polarised 
horn. 


Theoretical Background. 


[0152] When an object is illuminated by Electromagnetic 
(EM) waves the scattered waves in the time domain can be 
approximately represented by 


R(r)-ó(n) hr) 


[0153] where R is the scattered signal, I is the incident 
signal, h is the “impulse response” of the target, and 
represents the mathematical convolution operation, 

[0154] Ifthe incident signal is a sharp pulse approximating 
the delta function then the scattered (detected) signal 
approaches 


R(D)=8HDOO=h(0) 


[0155] So the scattered signal approximates the impulse 
response of the object. 

[0156] This is useful because the impulse response “car- 
ries” unique and aspect independent information about the 
electromagnetic response of the target, 


Phase (Aspect) 


Dependent) 
৩০ 
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Independent) Independent) 


[0157] Rewriting the late time response in exponential 
form we obtain 


A(t) = » Cm exp((—ay + 1278) + Cm CRD Cm — 2AVm)t 


m-l 


[0158] where the argument of the exponential terms are 
the complex natural resonances of the target—these are 
aspect independent. 

[0159] According to some embodiments of the invention, 
the Late Time Response (LTR) of the target 116 is recorded 
and Generalised Pencil Of Functions (GPOF) method is 
used to determine the complex natural resonances (poles)— 
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[0160] These resonances or poles are used to identify the 
presence of a particular (threat) object by comparing the 
poles of with those of measured objects and looking for a 
correlation. 
[0161] Requirements 
[0162] Amplitude and phase detection—not just power 
detection (e.g. use VNA) 
[0163] Ultra wideband (UWB) frequency to give suffi- 
cient time resolution to be able determine the frequency 
and decay times from the LTR 
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[0164] Where the microwave frequency is scanned from a 
lowest frequency v, through to a highest frequency v,, In 
currently preferred embodiments of the system, this is 
1/(17.5 GHz), i.e. 0.06 ns. 
[0165] Brief outline of technique: 
[0166] Detect scattered signals in both cross polarised 
and co-polarised orientation, 
[0167] Subtract stable background and transform to 
time domain (DIFFT) 


[0168] Time gate to isolate LTR and noise sample from 
signal 

[0169] DFFT LTR and noise sample data and subtract 

[0170] DIFFT to reconstruct “noise adjusted" LTR 

[0171] Use Generalised Pencil Of Function or other 


method to obtain complex natural resonances and apply 
statistical analysis of multiple acquisitions 


LTR Detection— Detailed Methodology 


[0172] Referring to FIG. 18A, a typical response of human 
body with concealed handgun is shown. The two traces 
1802, 1804 are for the two channels ofthe dual polarised TX 
horn 106' and RX horn 108' of FIG. 16. The early time 
response (BTR) and the late time (LTR) response are clearly 
distinguishable. The early time response is characterised by 
high frequency data while the low frequency oscillation after 
this is the late time response. The Late Time Response is 
aspect independent. In accordance with embodiments of the 
invention, the LTR must be separated from the ETR, as will 
be discussed further hereinafter. 


[0173] This is further illustrated in FIGS. 18B and C—the 
former is for a person without a handgun, and the latter for 
a person carrying a concealed handgun. 


[0174] FIGS. 19A to 19C show the process for detecting 
objects using LTRs, and in accordance with embodiments of 
the invention. 


[0175] As seen in FIG. 19A, in step s1902, an ultra 
wideband frequency sweep (at microwave frequencies) is 
taken, using the system of FIG. 16. Amplitude and phase 
data for scattered fields with (1) target and (ii) background 
(no target) are recorded in the frequency domain. 


[0176] In step 1902, the Late Time Response of a metallic 
object 1602 (FIG. 16) on a target 116 is excited. The object 
may be a benign object, or may be a threat object such as a 
handgun or knife, concealed on the human body 116. Due to 
the different LTR responses of the objects, it is possible to 
determine whether a threat object is present—this 1s illus- 
trated in FIG. 20C. The handgun is a more complex object 
and clearly displays a superposition of clamped resonances, 
compared with the simple objects of FIGS. 20A (spanner/ 
wrench) and 20B (Allen key). 

[0177] The excitement of LTRs is achieved by illuminat- 
ing the target 116 by a microwave source with frequency 
which is scanned from a lowest frequency v, through to a 
highest frequency v; to simulate illumination with a broad- 
band (frequency content) and hence short time length, 
electromagnetic pulse. v, is typically (but is not restricted to) 
0.5 GHz whereas v; is typically 18 GHz. The wavelength of 
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the lowest frequency must exceed the size of the object 1602 
in order to stimulate the fundamental mode (see further 
below). 


[0178] The microwave output has variable (user deter- 
mined) polarization state which is obtained via— 


[0179] (1a) Rotating a linear polarising filter 1610 imme- 
diately in front of the dual polarised microwave horn 106' 
used for target illumination (FIG. 16). The dual polarised 
horn 106' gives microwave output with orthogonal linear 
polarisation state; 


[0180] (1b) Introducing a time delay between the two 
input ports of the dual polarised microwave horn 106' used 
for target illumination in such a way as to give a time delay 


[0181] ‘Av between the orthogonal polarised components, 
where v is the frequency of microwave output; 


[0182] (2) Using a spiral free wave antenna (not shown) to 
give an elliptically polarised microwave output. 


[0183] The polarised microwave source with frequency v 
is scanned from a lowest frequency v,, through to a highest 
frequency v,. The frequency bounds of the system are 
selected so that the fundamental (lowest possible resonant 
frequency) and several higher order resonances of all threat 
objects 116 to be identified are encompassed within this 
range so that their excitation is possible. The natural reso- 
nant frequencies of objects 116 are related to the object’s 
largest linear dimension in as much as the fundamental 
frequency is—c/21, where | is the object's largest linear 
dimension and c is the velocity of light in free space. The 
upper bound is determined only by the apparatus used and 
should be as high as is practicably achievable since it is the 
bandwidth of the system that determines the resolution of 
time data obtained. 


[0184] The microwave receiving antennae 106', 108' (FIG. 
16) are connected to a data capture device(s) (detection 
electronics 1604, such as a Vector Network Analyser 
(VNA); FIG. 16) that is capable of measuring the complex 
amplitude (containing the magnitude and phase information 
of the scattered microwave field) of the scattered electro- 
magnetic return so that both amplitude and phase informa- 
tion can be utilised). The complex, frequency domain scat- 
tered signals of (in this case a human possibly carrying a 
concealed weapon) target and background (no target) are 
captured, Since the background is nominally constant this 
can be re-measured infrequently to ameliorate the effects of 
drift in the measuring equipment etc. 


[0185] Returning to FIG. 19A, subtraction (s1904) of the 
background response from the target response 15 carried out 
to give a complex, frequency domain signal which charac- 
terises the target only and is referred to as the target 
response. 


[0186] The target response data is then subjected (step 
s1906) to Inverse Fast Fourier Transform (IFFT) to give the 
time domain, scattered signal, with a time resolution At 
given by 


At= 





VH — VL 


[0187] The time domain signal so obtained is then time 
gated or ^windowed" (step s1908) to divide the time signal 
into four distinct sections (see FIG. 21) which correspond to 
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[0188] 1. The time period before the excitation 
[0189] 2. The time period during interaction with the 
propagating, exciting electromagnetic field 
[0190] 3. The time period immediately after excitation. 
[0191] 4. The time period after the LTR descends into 
the noise level 

[0192] FIG. 21 illustrates this separation, by showing the 
total time response for a Glock 17 handgun. The weapon is 
suspended in air. 
[0193] This separation of the time response into sections is 
achieved by either of two techniques—Method A and 
Method B, as discussed further below. 
[0194] Method A is illustrated in more detail in FIG. 19B. 
Here, sampling discrete, overlapping sections of the time 
response and applying a Fast Fourier Transform to each of 
these time samples. The frequency domain information thus 
obtained is Used to identify the beginning and end of the 
early time response as this section contains large amplitude, 
higher frequency data and the start of the LTR, as this 
section contains lower amplitude lower frequency data (can 
we show this with a figure of sample output?). Finally, the 
position at which the LTR descends into the random noise 
level of the signal is identified by correlation with the 
portion 1. 
[0195] First, at step s1930, a sliding time window is used 
to analyse (by FFT) spectral content of discrete windowed 
portions of time domain target response. The sliding win- 
dows have three user selected options—the segment size 
(time), the time shift applied to the window to “slide” it and 
whether the window is a step function or Gaussian in form, 
Typically both the window width and the time shift are <1 
ns and exact values are empirically determined. 
[0196] Next, at step s1932, both amplitude and frequency 
content of the spectral data (time domain target response) are 
analysed for each small time segment (window). This data is 
the used to determine where the Early Time Response starts 
and ends and where the Late Time Response starts and 
descends into noise (ends). 
[0197] As seen in step s1934, for each window, thresholds 
(Ta, T,) are applied to both amplitude (A) change and 
frequency (v) change compared to previous window. The 
ETR has large amplitude at high frequency whereas the LTR 
has large amplitude at lower frequency. The time periods 
before the ETR and after the LTR have low amplitude at all 
frequencies. The threshold for amplitude change and the 
threshold for frequency change are determined empirically 
by the user beforehand. 
[0198] Thus at step 51936, the position where Early Time 
Response starts is derived from (ampl. chge>T, and v 
chge>T,). 
[0199] Next, at step s1938, the position where Early Time 
Response ends/Late Time Response starts is derived from 
(ampl. chge<TA and v chge>T,). 
[0200] Finally, at step s1940, the position where Late Time 
Response descends into noise (ends) is derived from (ampl. 
chge»T , and v chge<T,). 
[0201] Method B is illustrated in more detail in FIG. 19C. 
This subprocess proceeds as follows. 
[0202] First (step s1942, the absolute maximum value in 
time domain target response is located, The absolute maxi- 
mum value is used to determine the position of the ETR 
(s1944) 
[0203] This Method, used two time delays (1 and 
2)—these may for example be retrieved from a database 
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(s1946). Time delay 1 (0) is chosen so that ETR from a 
typical human body is removed. This is done by applying the 
formula 


t=2D/c 


, where D is the width of a typical human body. Time delay 
2 (tj) is chosen so that a typical LTR response would have 
been attenuated into the noise level after this time, and is for 
example —5 ns. 


[0204] Next, the time domain target response is sampled 
(step s1948), starting from time delay 1 after time position 
of absolute maximum value and continuing for a time length 
equal to time delay 2. Finally, the sampled data is used 
and/or store as the LTR. 


[0205] Once the time gating (step s1908) has been per- 
formed, a pole extraction technique is applied to the derived 
LTR. 


[0206] Referring briefly to FIG. 21, the latter portion (3) is 
the Late Time Response. (LTR) and this complex data is 
defined as that immediately after the scattered response or 
Early Time Response (ETR). The LTR contains aspect 
independent information of conducting targets 116 carried 
on the body providing these conducting targets have natural 
resonance(s) which exist between the frequency bounds v,, 
and v,, with decay time(s) (lifetimes) that are sufficiently 
long that there is discernible amplitude above noise after the 
ETR of the body has passed. 


[0207] The LTR data S[n] i, expected to be of the form 


s[n] = GÈ Cm EXp(ZmnAt) + GÈ Cn 61002257451) + N [n] 


m-l m-l 


Where there exist M natural resonances Z,, between the 
frequencies v, and v,, there is assumed to be inherent noise 
in the system N. 


[0208] The natural resonances or poles Z,,=a,,+i2nv,, are 
aspect independent and are to be extracted in step s 1910 as 
precisely as is possible in order to confirm the presence or 
absence of a particular threat object 1602 whose natural 
resonances are known (either by measurement or numerical 
simulation) a-priori. The complex amplitudes C,, are highly 
aspect dependent and are not utilised explicitly in the 
determination of the presence or absence of a threat object. 


[0209] Depending on the embodiment, the poles are 
extracted using either the Generalised Pencil Of Functions 
method (Matrix Pencil method) or a Genetic Algorithm is 
implemented for this purpose. The complex poles and their 
associated complex amplitudes (residues) are thus extracted. 
The Generalised Pencil Of Functions method is well known 
to persons skilled in the art, and will not be discussed in 
detail here, see “Generalized Pencil-of-Function Method for 
Extracting Poles of an EM System from Its Transient 
Response" Yingbo Hua and Tapa K. Sarkar, IEEE Transac- 
tions on Antennas and Propagation, Vol. 37, No, 2, February 
1989. 


[0210] FIG. 22 illustrates poles extracted from LTR data 
using General Pencil of Functions method, varying model 
order from 1 to 10, for (a) a spanner, (b) an Allen key and 
(c), a handgun. It will be apparent that the position of the 
poles (circled) is different in each case. As the (approximate) 
position of poles for a number of threat objects is know, the 
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determination based on the poles in FIG. 22C enables the 
presence of a weapon to be determined, as discussed below. 
[0211] In an alternative embodiment, pole extraction (step 
$1910 in FIG. 19A) is performed using a Genetic Algorithm. 
For this case, and pseudocode for an exemplary algorithm is 
provided in Appendix C. 

[0212] A genetic algorithm and/or differential evolution- 
ary algorithm is a well established method of obtaining the 
parameters necessary to find the closest fit between the 
observed (experimental) Late Time Response (previously 
described as a sum of exponentially decaying sinusoidal 
functions or natural resonances) and the mathematical func- 
tion that describes them, (given earlier) [In-Sik Choi et. al. 
“Natural Frequency Extraction Using Late-Time Evolution- 
ary Programming-Based CLEAN”, IEEE Transactions on 
Antennas and Propagation, Vol. 51, No. 12, Dec. 2003]. In 
[In-Sik Choi et. al] the decaying natural resonances func- 
tions are fitted one at a time to the late time response and 
then subtracted from the original data. This makes the 
differential evolution algorithm task somewhat easier as it 
limits the number of parameters to four, amplitude, phase 
and the real and imaginary parts of each pole—a,,+i2nv,, per 
iteration. In a preferred embodiment, the maximum number 
of waveforms is fitted simultaneously, which is a more 
difficult task computationally as there are more permutations 
of possible solutions to explore. In a algorithm according to 
a preferred embodiment the crossover/mutation operator 
from differential evolution referred to in [Price and Storn 
1997] is used, followed by a tournament selection [Price and 
Storn 1997] to find the fittest chromosomes, although this 
does not preclude other least squares minimisation methods. 
It has been found that this approach is able to successfully 
fit the 20 parameters necessary to describe 5 damped sinu- 
soidal functions typically used when analysing responses 
from concealed guns and other weapons. 

[0213] In order to determine the presence of a weapon, 
two of the four parameters for each damped sinusoid (fre- 
quency find decay) are stored. Typically a handgun will 
require at least two damped sinusoids to classify it. If 
multiple data sets are acquired then a cluster of frequencies 
and damping factors can be used to identify the weapon (see 
FIG. 22C). In order to classify the weapon at least two 
frequency/damping factor combinations can be used. 
Because the extracted parameters will never be exactly 
reproducible due to noise and other experimental artifacts, a 
clustering technique such as fuzzy c or k means clustering 
can be employed to find the centre of each cluster, or a self 
organising map type neural network will distribute neurons 
such that a particular set of neuron will be activated if a 
particular cluster pattern is presented to the network. 
[0214] Returning to FIG. 19A, following pole extraction 
(s1910), the poles are filtered (s1912). This means that poles 
are discarded if the magnitude of their amplitudes are 
smaller than a user set threshold or if the damping constants 
(real part of pole) are positive or if the frequency (imaginary 
part of pole) is negative, A user set maximum and minimum 
is also applied 

to both the damping and frequency parts of the pole, If the 
pole lies outside of this space then that pole is ignored. This 
discrimination space is selected from experimentally mea- 
sured poles of likely threat objects, i.e. by measuring the 
poles of a selection of handguns on a body. 

[0215] Steps 1902 to 11912 are repeated, as shown by loop 
s1914 a set number of times— decided by user set parameter. 
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The number of loops is determined by the likely time a 
person being interrogated would spends within the active 
area (e.g. see corridor of FIG. 17) of the system, and this will 
depend on individual system configurations, 

[0216] Next, following pole filtering (step s1912), the 
poles are stored S1916). 

[0217] Then, at step $1918, the pole data is compared with 
library (1920) of measured poles for targets of interest, 
stored in a database. 

[0218] The comparison is carried out by finding the closest 
matching library poles to those measured and then comput- 
ing a root mean square error where the damping space and 
frequency space are weighted by experimentally determined 
values. This weighting is necessary as the pole position in 
damping space is more spread than in frequency space. The 
RMSE value is then compared to empirically determined 
threshold values to give a threat level based on the closeness 
of pole match. 

[0219] Finally, it a threat level decision (0-1) is obtained 
and a possible threat object/class determined (s1922). 
[0220] Embodiments of the invention include the follow- 
ing novel aspects and consequent advantages, 

[0221] 1. Ultra Wide Band (UWB) to give large fre- 
quency coverage and thus excite maximal number of 
resonances and to give high time resolution (short time 
span) in LTR 

[0222] sufficient that rapidly oscillating and quickly 
decaying resonances can be captured. 

[0223] 2. Robust auto separation of the early and late 
time domains—giving the LTR for pole extraction, 
partly facilitated by the use of Ultra Wide Band exci- 
tation to make the boundary more obvious. 

[0224] 3. Anechoic portal design to give low noise 
response data. 

[0225] 4. Multiple pairs of transmitter/receiver anten- 
nae for all round target interrogation. 

[0226] 5. Continuous or stepped, scanned polarisation 
state through mechanical or phase generated elliptical 
polarised output to optimise possibilities of coupling 
into aspect independent modes. 

[0227] (i) A type of Genetic Algorithm known as an 
Evolutionary program for processing multiple sweeps 
of LTR data and simultaneously extracting the complex 
natural resonance poles for target identification, form- 
ing clusters of poles. The use of Evolutionary program- 
ming does not preclude the use of more conventional 
techniques such as Pencil of Function methods. 

[0228] Gi) Cross polarised transmission and receiving 
antennae to give enhance discrimination between body 
ETR and LTR. 

[0229] (ii) The pre-excitation time domain data can be 
used as a measure of the effectiveness of the back- 
ground subtraction, since if there are no scattering 
surfaces present between transmitter and target this 
data should be of very small amplitude relative to the 
other portions. The Fast Fourier Transform of I the 
pre-excitation time portion can be used to improve the 
LTR data, by subtraction in the frequency domain and 
subsequent time domain reconstruction. 

[0230] (iv) Undesirable effects of non-uniform trans- 
mitting and receiving antenna response can be miti- 
gated by division of the frequency domain target 
response by the absolute value of the measured horn 
response in the frequency domain. 
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[0231] (v) Multiple transmitter/receiver pairs to give all 
round target coverage. 


APPENDICES 
[0232] 


APPENDIX A 


Procedure for collection and automatic analysis of threat object sensor 
Receive activation signal (e.g. operator activates scan button or subject 
triggers scan] 
For a predetermined number of sweeps do 

While full frequency range not scanned do 

Illuminate subject with radiation 

Step over frequency range 

Receive reflected radiation signals 

Perform transformation on received radiation signals to 

produce time domain or optical depth domain trace 
Store in a sweep channel 

End While 

Increment sweep channel 
End do 
Normalise time domain trace according to range 
Use Complex Fourier Transform (VNA mode) or Direct Fourier 
Transform (reflectometry mode) to convert to x-dimension to 
determine position of trace peaks. 











From x-positions use conversion factors L = 


[4 


JAF to determine corrected x-axis (optical depth). 


++++++++++++++++++717 





Perform transformation on received radiation signals by all 
the steps below: 


APPENDIX B.1 





[pseudocode for transformation for Swept reflectrometry] 
If (Technique 1 (Swept reflectrometry)) then 
Use Direct Fourier Transformed signals (fft of IER?) 
For each sweep channel do 
Set Lower and Upper bands LI L2 for useful optical depths 
(e.g. 10mm to 150mm depending on weapon size and orientation). 
Set Threshold for useful signal level above previously collected 
values for body alone. 
From LI to L2 do 
Store Signal above Threshold separately 
in vectors in arrayl 
Integrate Signals above threshold 
End do 
End do 
For each sweep channel do 
Correlate adjacent vectors and produce outputl 
Sum with previous outputl's 
Sum integrated signals above threshold 
End do 
Outputl is sum of all correlations between vectors in arrayl 
Output2 is sum of integrated signals above threshold for each sweep 
channel 
Outputl will be different for gun when optical depths change as subject 
moves in beam than for block of explosive stimulant which is of a 
similar thickness from different aspects. 
Outputl and Output2 are taken to Neural Network input (see FIG. 15). 
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Else If (Technique2 (Barrel tone detection)) then 
[pseudocode for transformation for Barrel tone detection] 
Use Direct Untransformed signals 17371 
For each sweep channel do 
For set of weapon calibers do 
For set of barrel lengths do 
Calculate onset (f0) for caliber 
Calculate chirped response for Length L 
Correlate ideal response 
cos(2n(1?- fo?) ^ (2L/c)4(q) with data. 
Store Correlation value 
End do 
End do 
Find best (lowest) correlation value and store in Output3 
End do 
Keep 10 and L and display 
Output3 goes to Neural network (see FIG. 15). 
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Else if (Technigue3 (Cross-polarization detection)) then 
[pseudocode for transformation for Cross-polariasation detection] 
Use Complex Fourier Transform signals (Eg) which give range 
information, from normal and cross polarized detectors 
Select Distancel which is first significant reflection above a threshold 
(or which is given by an independent range finding sensor) 
For each sweep channel do 
Apply a distance window of a given number of millimeters 
determined by database of responses from weapons. 
Select trailing edge of response (distance2) by adjusting 
distance window to when response falls below threshold 
Integrate response from non-pol detector within 
window 
Integrate response from cross-pol detector within 
window 
Sum to previous integrations 
Correlate response from non-pol detector with 
cross-pol 
Store and sum correlation 
End do 
Outputs 4 and 5 are sum of correlations and sum of integrations 
Outputs are taken to Neural Network inputs (see FIG. 15). 
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Else if (Technique4 (Late time response detection and resonant frequency detection)) then 


[pseudocode for transformation for Late time response] 
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Use Complex Fourier Transform signals (Eg) which give range information, from normal and 


cross polarized detectors 
For each sweep channel do 
Select Distance2 which is the output of Technique3 


Apply a distance window of a given number of millimeters determined by 


database of late time responses from weapons. 


Find the smoothed responses of the late time responses (see FIG 14) in normal 


and cross-pol detectors. 


Ascertain the exponential decay rate by a non-linear least squares fitting 


procedure. Output is stored in a vector 
End do 
Normalise outputs from the vector and form Output 6 


Use Complex Fourier transformed signals (Eg), from normal and cross polarized 


detectors 
For each sweep channel do 
For the entire sweep data 
If response level above a normalised threshold then 





Apply series of non-linear filters (e.g. MUSIC filter) with filter characteristics 
taken from a data base to look for a particular resonance Store the magnitude 


the resonance Apply a peak detection algorithm (E.G. zero crossing) Store peak 


locations 


For each sweep channel do 





Compare peak locations with known natural resonances for object 


Sum the differences between peak locations and natural resonances from data base. 


Output7 is sum of differences - to Neural Network. 


Keep the peak locations for display, which can indicate weapon type. 
End do 

End if // End of transformation techniques phase 
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Task for Differential Evolutionary program 
To minimise the sum of sguares error (SSE) between the observed 
and calculated late time response (LTR). 


LTR generated from 


1) < le 
Equation, s[n] = (5) > Cmexp(Z,nAt) + (5)১ Cm given in text 


m-l m-l 


For this algorithm POLE indicates frequency, phase shift, decay 
rate and amplitude. 

To choose the optimum number of poles by monitoring decrease 
in SSE after each increase in the number of poles 

To obtain the optimum LTR data by a SIMULTANEOUS fit of 
multiple poles unlike the technique of Choi 'et al who find 
optimum parameters for a SINGLE pole and then extract this 
rom the original data set ITERATIVELY until the desired 


number of poles is obtained. No criteria are given for the 








ermination of the iteration. 
ALGORITHM 
LOOP L from MINPOLES to MAXPOLES* 

LOOP K from 1 to MAXIMUM NUMBER OF GENERATIONS 





Create initial random population of solutions 
LOOP M from 1 to SIZE OF RANDOM POPULATION 
Calculate observed LTR with L poles according to M^ 








population data 
Perform SSE error calculation 
Store Mth SSE error 

END of loop M 





APPENDIX C-continued 





Store minimum SSE error and associated gene 
Mutate gene pool according to Ref, [2] 
END OF LOOP K 
END OF LOOP L 





*Note: 

MINPOLES typically 2 and MAXPOLES typically 5 

NEGLECT POLES WITH VERY SMALL AMPLITUDE 

NEGLECT POLES WHEN SSE STILL IMPROVING ACCORDING TO A USER 
DEFINED TOLERANCE 
STORE OPTIMUM FREQUENCY AND DECAY RATE 

Don't store AMPLITUDE AND PHASE - not used for eventual weapons classification. 
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1-93. (canceled) 


94. A system for remote detection of one or more dimen- 
sions of a metallic and/or dielectric object, comprising: 


a transmission apparatus, including a directional trans- 
mission element, configured to be pointed towards a 
target by a user, the transmission apparatus being 
configured to direct microwave and/or millimeter wave 
radiation with a first polarization in a predetermined 
direction, 


a non-imaging detection apparatus configured to detect 
radiation from an entity resulting from the transmitted 
radiation, the detection apparatus being configured to 
detect radiation which has the first polarization and to 
detect radiation which has a second polarization which 
is transverse to the first polarization, the detection 
apparatus being further configured to generate one or 
more detection signals in the frequency domain, and 


US 2019/0391256 Al 


a controller, the controller being operable to 
(i) cause the transmitted radiation to be swept over a 
predetermined range of frequencies, 

(11) perform a transform operation on the detection signal(s) 
to generate one or more transformed signals in the time 
domain, and 

(111) determine, from one or more features of the transformed 
signal, one or more dimensions of a metallic or dielectric 
object upon which the transmitted radiation is incident. 

95. The system of claim 94, wherein the controller is 
operable to monitor for anti-correlation between detected 
radiation which has the first polarization and detected radia- 
tion which has the second polarization. 

96. The system of claim 94, wherein the detected returned 
radiation which has the first polarization is analyzed in the 
time domain, and the detected returned radiation which has 
the second polarization is analyzed in the time domain. 

97. The system of claim 94, wherein the first polarization 
is detected using a first receiver horn and the second 
polarization is detected using a second receiver horn. 

98. The system of claim 94, wherein the first and second 
polarizations are measured using receivers which are ori- 
ented parallel to the direction of transmission of the radia- 
tion. 

99. The system of claim 94, wherein the first and second 
polarizations are measured using receivers which are ori- 
ented at right angles to the direction of transmission of the 
radiation. 
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100. The system of claim 94, wherein the controller is 
operable to initiate step (1) upon receiving an activation 
signal, the activation signal corresponding to a user input 
and/or detection of the presence of the entity. 

101. The system of claim 94, wherein the controller is 
configured to normalize according to a range one or more 
transformed signals. 

102. A system for remote identification of the caliber 
and/or barrel length of a metallic threat object, comprising: 

a transmission apparatus, including a transmission ele- 
ment, configured to direct microwave and/or millimeter 
wave radiation in a predetermined direction, 

a detection apparatus configured to receive radiation from 
an entity resulting from the transmitted radiation and to 
store an amplitude of the received radiation as a func- 
tion of frequency, and 

control circuitry, the control circuitry being operable to 

(i) cause the transmitted radiation to be swept over a 
predetermined range of frequencies, 

(11) perform a transform operation on the detection signal 
(s) to generate one or more transformed signals in the 
time domain, 

(iii) extract, from one or more features of the transformed 
signal, the onset of the oscillatory response; 

(iv) identify the barrel length of the metallic threat object 
from the frequency of onset of the oscillatory response. 
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